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ABSTRACT 

We present the first high time-resolution simultaneous X-ray, ultraviolet, and optical observations of 
X-ray bursts in UY Vol, the optical counterpart of the low mass X-ray binary EXO 0748-676, obtained 
with RXTE, HST, and Gemini-S. Strong reprocessed signals are present in the ultraviolet (a factor 
of 4) and optical (a factor of 2.5). These signals are lagged with respect to the X-rays and appear 
significantly smeared as well. The addition of far- ultraviolet coverage for one burst, together with 
the high quality of the dataset, allow much tighter constraints upon the temperature and geometry 
of the reprocessing region than previously possible. A single-zone black body reprocessing model for 
this burst suggests a rise in temperatures during the burst from 18 000 to 35 000 K and an emitting 
area comparable to that expected for the disk and/or irradiated companion star. The lags, a mean of 
4.0 s and range of 2.5 s, are consistent with those expected within the binary. The single-zone black 
body model cannot reproduce the ratio of optical to ultraviolet flux during the burst, however. The 
discrepancy, corresponding to underpredicting the optical by more than a factor of two, seems too large 
to explain with deviations from a local black body spectrum and more likely indicates that a range 
of reprocessing temperatures are required, as would be expected, with cooler regions not contributing 
to the UV. Comparable results are derived from other bursts, and in particular the lag and smearing 
both appear shorter when the companion star is on the near side of the disk as predicted. The burst 
observed by HST also yielded a spectrum of the reprocessed light. It is dominated by continuum, 
with a spectral shape consistent with the temperatures derived from lightcurve modeling. Some line 
enhancements are also seen, most prominently in Cm 1175 A. Taken as a whole, our observations 
confirm the standard paradigm of prompt reprocessing distributed across the disk and companion 
star, with the response dominated by a thermalized continuum rather than by emission lines. 

Subject headings: accretion, accretion disks — binaries: close - stars: individual: UY Vol 



1. INTRODUCTION 

The low-mass X-ray binary (LMXB) EXO 0748-676 
was discovered in 1985 as a transient X-ray source 
l)Parmar et al]ll98fll) and ra pidly associated wi th an op- 
tical counterpart, UY Vol ( Wa de et alJ 11985). Unlike 
most X-ray transients, however, it did not decay back to 
a quiescent state, but remained active and is now con- 
sidered part of the persistent LMXB population. This 
makes it an intriguing object for study as it appears to 
exist near the edge of stability, flipping between phases 
of quasi-stable activity and quiescence. The likely key 
to its behavior is that the system is currently held in a 
meta-stable high state by i) stabilization of the accretion 
disk against thermal instability by X-ray irradiation and 
ii) enhancement of the mass transfer rate from the com- 
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panion. UY Vol is thus an important object for under- 
standing the impact of X-ray irradiation upon the disk 
and companion star in LMXBs. 

UY Vol also has an inclination that is just right to 
give us a range of diagnostic tools to probe the accretion 
flow. The inclination is high enough that X-ray eclipses 
are seen re curring on a 3.82 h r period, and X-ray dips 
also occur ijParmar et al.lll986]) . The eclipses are sharp, 
indicating that the neutron star itself is being eclipsed, 
and hence that it is visible outside of eclipse. UY Vol is 
therefore not an accretion disk corona (ADC) source; an 
ADC is likely present, but is not the dominant source of 
observed X-rays. 

To better understand the accretion structure and effect 
of irradiation in this LMXB, we have performed a multi- 
wavelength study using HST for UV rapid spectroscopy, 
RXTE for X-ray data, Gemini-S for rapid optical pho- 
tometry, and the Cerro Tololo 4 m Blanco telescope for 
optical spectroscopy. We discuss here the analysis of sev- 
eral X-ray bursts seen during the simultaneous coverage. 
Subsequent papers will address the UV and optical emis- 
sion line spectra ijPearson et alJ 120061 hereafter Paper 
II), and the spectral energy distribution, and multiwave- 
length orbital light curves (Paper III). 

Type I X-ray bursts, seen in low-mass X-ray binaries 
(LMXBs) with neutron star primaries, are due to explo- 
sive thermonuclear burning of accreted material on the 
surface of the neutron star. They involve a large increase 
in the X-ray flux, by a factor often or more, on timescales 
of a few seconds. As well as providing insights into the 
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TABLE 1 

Log of UV, optical, and X-ray high time resolution 

OBSERVATIONS. 



Facility Instrumentation Start date UT range Total 

time (s) 



HST 


STIS, G140L 


2003 Feb 18- 


19 


20:06- 


-00:06 


13470 


HST 


STIS, G230L 


2003 Feb 19 




00:17- 


-00:31 


800 


HST 


STIS, G140L 


2003 Feb 19 




00:53- 


-04:52 


13470 


HST 


STIS, G230L 


2003 Feb 19 




05:04- 


-05:17 


800 


Gemini-S AcqCam, V 


2003 Feb 18 




04:19- 


-05:32 


9600 


Gemini-S AcqCam, V 


2003 Feb 19 




01:08- 


-08:38 


18960 


RXTE 


PCA 


2003 Feb 14 




00:08- 


-09:32 


10300 


RXTE 


PCA 


2003 Feb 15 




01:27- 


-09:12 


8238 


RXTE 


PCA 


2003 Feb 17 




17:21- 


-18:23 


3520 


RXTE 


PCA 


2003 Feb 17- 


IS 


20:30- 


-08:13 


11751 


RXTE 


PCA 


2003 Feb 18 




12:17- 


-13:19 


1792 


RXTE 


PCA 


2003 Feb 18- 


19 


20:11- 


-09:23 


13706 


RXTE 


PCA 


2003 Feb 19 




15:07- 


-16:09 


3616 



conditions on the surface of the neutron star, the sud- 
den flash lights up the whole binary system, and signif- 
icant reprocessed bursts, with an amp litude of a fa ctor 
of a few, are seen in the optical (see iHvnesI 120051 and 
references therein). Since the reprocessed bursts are a 
large amplitude signal rather than a small perturbation, 
the non-linearity of the optical response (due to observ- 
ing band-limited rather than bolometric fluxes) is clear. 
Thus lightcurves in different bandpasses exhibit signif- 
icant differences, breaking the degeneracy between the 
temperature and emitting area of the reprocessing region 
that exists when only s mall perturbations are considered 
ijLawrence et al.l|l983j) . This makes the bursts a very 
powerful tool for ap plying echo mapping techni ques in 
an X-ray binary (e.g. lO'Brien et al.ll2002[lHvnesll2005|) . 

In spite of this potential, the unpredictable nature of 
the bursts means that relatively few simultaneous obser- 
vations exist of any X-ray bursting source, many twenty 
years old, and no UV observations of bursts have been 
obtained. We present here the first simultaneous X-ray, 
UV, and optical observations of a burst in UY Vol, in- 
cluding fast far-UV spectroscopy. We also present several 
more bursts observed with only X-ray and optical cover- 
age. 

2. OBSERVATIONS 
2.1. HST 

HST observations used the Space Tel escope Imaging 
Spectrograph fSTIS: iProfitt et all 120021. with the far- 
UV MAMA detector and the G140L grating, and are 
detailed in Table □ TIMETAG mode yielded a stream 
of detected events, with 125 /zs precision, which could be 
used to reconstruct spectra for any desired time-interval, 
as well as high time-resolution lightcurves. The obser- 
vations were timed such that the target was within the 
continuous viewing zone (CVZ). Consequently we were 
able to observe over about 9 hrs with only small gaps for 
wavelength calibrations and mode changes. This covered 
two complete binary orbits. 

To obtain a lightcurve, we extracted source counts from 
a 50 pixel wide window centered on the target and back- 
ground from two similar regions on either side. Geo- 



coronal Hi Lya and Ol 1304 A lines were masked out of 
both source and background regions, as were the extreme 
ends of the detector. It was necessary to use this back- 
ground subtraction procedure as the background was 
found to be larger than the nominal global 7cnts~ 1 dark 
current and time dep endent. This inter mittent 'glow' 
is not easily modeled l)Prontt et alJ 12002]) so is best re- 
moved empirically. We did this by subtracting a poly- 
nomial fit to the background lightcurve. Source counts 
out of burst were ~ 30 — 100 s" 1 and the estimated 
background within the source window was < 7s _1 ; for 
comparison the nominal dark current should only be 
~ 0.3 s -1 within this window. 

We used 1 s time-resolution to approximately match 
the Gemini-S opt ical photometry. As discussed by 
iHvnes et af]l|2003D . #ST/STIS absolute timing accuracy 
is uncertain at a level of up to a few seconds, so we did 
not attempt precise barycenter corrections and allowed 
the zero point of the HST timing to be a free parameter 
in subsequent analysis. 

Flux calibration of the lightcurves was derived from 
the average spectra, and takes advantage of the redden- 
ing measurement (E(B — V) = 0.06 ± 0.03 9 ) possible 
with complementary near-UV data (see Paper III). For 
each G140L spectrum we compared the observed counts 
as a function of w avelength with the dereddened fluxes 
l)Fitzpatricld ^99) to determine the effective sensitiv- 
ity, hence defining our far-UV bandpass. We then used 
this sensitivity function to derive weighted average wave- 
lengths and fluxes. These averages were consistent to 
< 1 A in wavelength, and < 1 % in flux between the six 
G140L spectra. The effective wavelength derived (for our 
spectral shape) was 1388 A. The uncertainty in the cal- 
ibration of observed fluxes is 4 % (from the documented 
flux calibration of low-resolution STIS/MAMA modes; 
IProfitt et al]l200*2T) . However, the dominant term is due 
to the uncertainty in dereddening - this introduces a 
±25 percent uncertainty in the absolute calibration of the 
far-UV lightcurves, and ±16 percent in the relative cali- 
bration of the UV and optical fluxes (since these errors 
are corr elated). If the ext inction curve differs from the 
assumed iFitzpatrickl l)1999f) form then the error could be 
larger. 

2.2. Gemini-S 

On the nights of February 17-18 and 18-19, we used 
the Acquisition Camera (AcqCam) on Gemini-South to 
obtain fast V band optical photometry (see Table for 
details). Conditions were photometric, with realized im- 
age quality on target mostly ~ 0.9arcsec, although this 
degraded to ~ 1.2 arcsec in the latter part of the second 
night. For practical data acquisition it was necessary to 
break each night into series of (usually) 4000 images at 
a time. The series obtained, totaling 37000 images, are 
listed in Table [2 

One of the great strengths of this instrument for fast 
photometry is the fast readout. We used the camera 
windowed and binned (2x2 giving 256 x 256 0.24 arcsec 

9 Note t hat this is s ubstantial ly below the previously quoted 
value of 0.4 iSchocnibs & Zocschingcr 1990; Liu et al. 2001). Since 
that value was based on shifting the object to lie on the main- 
sequence in a color-color diagram, it is not expected to be reliable 
for an LMXB 
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binned pixels). Other data acquisition modifications 
were made to further minimize the dead-time between 
images. This dead-time was not absolutely constant, but 
was usually 0.305-0.310 s, with occasional (less than one 
per few hundred) glitches to as high as 0.44 s. The precise 
time-stamps were taken from when the image completed 
writing to disk, so the start time of the subsequent image 
(relative to the first in a series) was known. The start 
time of the first image is taken from a GPS based clock. 

Basic reductions including bias removal, subtraction of 
the significant dark current, and flat-fielding were done 
using the AGREDUCE script within iraf 10 . The win- 
dowed field includes one star much brighter than UY Vol, 
which provided a high fidelity reference star and two 
other comparisons (one brighter, one fainter) which were 
used to verify the accuracy of the results. We experi- 
mented with both small aperture photometry using iraf 
and optimal photome try as implemented in the Starlink 
PHOTOM package Csee lNavlorlll998l for the algorithm this 
is based on). We found negligible difference between 
the two methods provided the aperture (for unweighted 
photometry) was optimized. We opted to use the re- 
sults from the optimal algorithm as this explicitly ad- 
justs the weighting per image and so should be more 
robust against changes in seeing. We used the compar- 
ison stars to verify that the standard deviations of the 
resulting lightcurves are dominated by the formal errors 
in the photometry. The formal errors for UY Vol were 
typically ~ 2.5 percent per 0.8 sec exposure. 

Absolute calibration was done with respect to two 
stars from SA104 (Landolt 1992) on the first night. 
We derive magnitudes of V = 13.63 for the refer- 
ence star and 17.4-18.0 for UY Vol (excluding bursts). 
The magnitude range derived for UY Vol is within 
the range observed by earlier studies (e.g. 17.1-18.1; 
Ivan Paradiis. van der Kris, fc Pederserj fl988). The or- 
bital lightcurves will be considered in a subsequent pa- 
per. 

Finally, we corrected the observed magnitudes for 
interstellar extinction (assuming E(B — V) = 0.06 
based on the 2175 A feature in our HST data) and 
converted to fluxes with the conversi on constant of 
lFukugita~ Shimasa ku. fc Ichikawal ©95). 

2.3. RXTE 

X-ray observations were obtained with the Rossi X-ray 
Timing Explorer (RXTE) in two blocks timed to coincide 
with optical spectroscopy (Paper II) , and the optical and 
UV observations described above. A log is presented in 
Table H] 

X-ray lightcurves of the bursts were recovered from 
a 64 channel even t mode of the Prop ortional Counter 
Array (PCA). See lJahoda et all l)2006|) for discussion of 
the current status of the PCA. Since the source counts 
are high during a burst, and the timescale is short, we 
used all PCUs which were switched on, including PCU0. 
We initially processed lightcurves from each PCU sep- 
arately, however, to check that no background flares 
were present in PCU0, before combining them. Back- 

10 IRAF is distributed by the National Optical Astronomy Ob- 
servatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 



ground lightcurves (with 16 s time-resolution) were con- 
structed from the L7/240 faint source combined mod- 
els dated 2002 Feb 1. This should allow the most reli- 
able background subtraction from the pre-and post-burst 
lightcurves. Since the L7 rate is modified for bright 
sources, however, this model lightcurve is contaminated 
during the burst. We therefore interpolated between the 
pre- and post-burst modeled background rates to define 
the burst background and subtracted this. We note that 
in practice the background subtraction has little impact 
on the results described here as we perform an addi- 
tional empirical subtraction of the persistent (pre- and 
post-burst) flux. Hence more refined background models 
would not significantly affect our results. 

Lightcurves were initially extracted for 2-5 keV, 5- 
12keV and 12-60keV bandpasses. The burst profiles 
varied significantly with energy, so it is important to be 
careful in choosing an appropriate lightcurve for a de- 
convolution. We therefore also constructed a 2-20 keV 
integrated flux lightcurve. To do this, we estimated ap- 
proximate per-channel conversions from count rates to 
fluxes using XSPEC and appropriate response matrices. 
We used a black body fit to the average burst spectrum 
to define these conversions to ensure that the relative 
weighting within each channel was approximately cor- 
rect. We then applied these conversions to each channel 
and summed the fluxes. This will be somewhat noisier 
than a straight sum of count rates, but represents our 
best estimate of the evolution of the flux irradiating the 
disk and the companion star. Based on the black body 
fits to the spectra, the 2-20 keV bandpass accounts for 
most of the burst flux (> 90 %), so additional bolometric 
corrections are ignored. 

3. SYSTEM PARAMETERS 

Any attempt to quantitatively model observations of 
UY Vol requires system parameters. UY Vol does not 
yet benefit from dynamical estimates, so we must use 
more indirect constraints. 

Fortunately it is eclipsing, and precisely defined X-ray 
eclipses provide exquisitely detailed measurements of the 
orbital period an d duration of the neutron star eclipse 
( Wolff et al.l2002j) . Through these eclipses, we know that 
the inclination must be high, and the relationship be- 
tween mass ratio and inclination is well defined. Assum- 
ing recent eclipse durations of 497.5 ± 6.0 s, we obtain the 
relationship shown in Fig. ^ We can attempt to further 
constrain the available parameter space in other ways. 
The fact that we see the sharp neutron star eclipse at 
all indicates that UY Vol is not an accretion disk corona 
(ADC) source. This means the inclination cannot be too 
high as our line-of-sight must pass over the disk rim. The 
disk rim height is not known directly, but must be greater 
than that expected from hydrostatic support alone. This 
is a very weak constraint, however; for example, for a 
disk half-thickness of 0.03, we only require i < 88.3°. 
This constraint is shown in Fig. ^ 

Other constraints are more model-dependent. We can 
make plausible estimates of the range of mass ratios 
likely, although values outside of this range are still pos- 
sible. The companion star to UY Vol is probably above 
the hydrogen burning l imit and non-degenerate. This 
implies M 2 > 0.07 M llChabrier fc Baraffd 120001). As- 
suming a 1.35 M neutron star (jThorsetTTT Chakrabartvl 
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Model 2: 3=0.20, 1=78 
MB1: Phase=0.20 



Model 2: (3=0.20, 1=78 
MB2 i Phase=0.38 




Fig. 1. — Relation between inclination and mass ratio based on 
the widths of X-ray eclipses. 



1999) then yields a minimum mass ratio of 0.05. This 
limit is soft, as a more massive neutron star is certainly 
possible if it has accreted significant mass from the com- 
panion (as seems likely in 2S 0921-630; Shahba z et al.l 
l200ll,Tonker et al.l2005fl . At the high mass-ratio end, we 
can assume that the companion star is underdense com- 
pared to a main-sequence star that would fill the Roche 
lobe. Again assuming a 1.35 M© neutron star this yields 
the upper limit on the mass ratio shown. This could be 
increased a little for a somewhat lower mass neutron star, 
and further if the companion is actually overdense. This 
is possible if the companion had undergone nuclear evo- 
lution, before losing most of its envelope in mass transfer 
ijSchenker fc Kind2002llHaswell et all2002|) . In this case 
the neutron star might also be expected to be more mas- 
sive. Consequently, there is probably not much scope for 
a mass ratio higher than 0.4 in this system. 

It should be noted that the constraints considered so 
far are effectively the same as obtained bv lParmar et alJ 
(1986), namely 75° < i < 82°. As discussed above, the 
limits are rather soft, and a slightly larger range is pos- 
sible with a companion star and/or neutron star with 
extreme properties. 

To obtain other parameters of interest, we synthesize 
a population of possible binaries. We assume the known 
orbital period, and the relationship between mass ratio 
and inclination given in Fig. 2] We consider neutron star 
masses with an asymmetric G aussian distribution Mi = 
1-351qq 4 M©, i.e., following l)Thorsett & Ch akrabartv 
1999), but allowing for a higher mass tail due to mass 
transfer. We assume a uniform distribution of compan- 
ion star masses between 0.07 M©, and the mass at which 
main-sequence density is reached. 

We derive a relatively narrow range of binary sepa- 
rations, (1.03 ± 0.05) x 10 n cm (at 90% confidence), 
since the binary period is known and there is not a 
large uncertainty in the total system mass. Disk pa- 
rameters are more uncertain; the tidal truncation radius 
is (0.50 ± 0.05) x 10 11 cm and the projected area of a 
flat disk would be (1.5 ± 0.2) x 10 21 cm 2 . The projected 
area of the companion star, with a spherical approxima- 
tion, would be (2.0 ± 0.9) x 10 21 cm 2 . While in general, 



Model 2: q=0.20, i=78 
MB3: Phase=0.48-0.51 




Fig. 2. — System geometry at the phases where bursts were 
observed for model 2. 



the companion is expected to subtend a larger projected 
area than the disk, only a phase-dependent fraction of 
this will be X-ray heated, so this should be considered 
an upper limit on the area of the luminous regions of the 
companion, which could be significantly less than that of 
the disk. 

We can estimate light travel times in the same way. 
The maximum lag from the pole of the companion star 
would be (6.8 ± 0.3) s at phase 0.5, although most of the 
heated inner face would of course be a little shorter than 
this. Lags from the disk would be expected to span a 
range from zero to (3.3 ± 0.3) s assuming it extends to 
the tidal truncation radius. 

We define several test cases consistent with Fig. ^ 
model 1 has q = 0.08, i = 82°. Model 2 has q = 0.2, 
i = 78°, and model 3 has q = 0.34, i = 75.5°. Fig. El 
shows a schematic view corresponding to model 2 at 
phases when bursts were observed. 

An important factor in considering irradiation of the 
companion star is the opening angle of the disk, 0. For 
our purposes, this defines the height of X-ray absorbing 
material, which may be above the optical photosphere, 
and may not be in hydrostatic equilibrium (for example 
material thrown up from the stream impact point or local 
flares). Values deri ved from other objects have generally 
been rather high. Ide Jong, van Paradiis. & Aug usteiinl 
( 1996) derived (3 = 12° and cite other authors who ob- 
tained a range of 6-14°. In our case, the highest val- 
ues are ruled out in some models, as UY Vol is not an 
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TABLE 2 

Bursts detected during our observations 

1.0 



Burst ID Wavelengths UT date and time Phase 



XB1 


X-ray 


2003 Feb 15, 05:27:06 0.22 


XB2 


X-ray 


2003 Feb 18, 01:15:16 0.95 


MB1 


X-ray, opt. 


2003 Feb 18, 06:01:38 0.20 


MB2 


X-ray, opt., UV 2003 Feb 19, 01:49:12 0.38 


MB3a 


X-ray, opt. 


2003 Feb 19, 06:01:36 0.48 


MB3b 


X-ray, opt. 


2003 Feb 19, 06:09:08 0.51 



ADC source - we see neutron star eclipses and promi- 
nent bursts, hence we do observe the neutron star di- 
rectly. Thus the opening angle must be less than 90 — i, 

1. e. less than 12° for model 1 and less than 8° for model 

2. This is not a constraint for model 3. Note that the 
presence of eclipses also indicates that at least some of 
the companion star must be exposed to direct radiation 
from the neutron star; it cannot be fully shielded by the 
disk. 

A final important parameter is the distance to UY Vol. 
The best indicator of this is in neutron star LMXBs is 
the peak flux observed during radius expansion X-ray 
bursts as th is is believed to be an approximate stan- 
dard candle ijKuulkers et al.l 2003*1. Three radius ex pan- 
sion bursts were reported by iGottwald et all l|1986l ) and 
IJonker fc Nelemana 1120041) used the se, together with the 
calibration of l|Kuulkers et alJl20Q^l . to derive a distance 
range to UY Vol of 6.8-9.1 kpc, with the low end of the 
range corresponding to burning of material of normal 
composition and the upper end t o hyd rogen poor ma- 
terial. More recently, iWolff et aD 1)20051) used a brighter 
radius expansion burst seen by RX TE to obtain distances 
of 5.9±0.9kpc or 7.7±0.9 kpc for hydrogen rich and poor 
bursts respectively. The latter authors argue that the 
variation in brightness likely reflects variable obscuration 
and hence that their brighter burst gives a more reliable 
distance measurement. We will thus use the midpoint 
of the latter estimates for this work, 6.8 kpc. Allowing 
for the composition ambiguity, all luminosities and esti- 
mated areas quoted in this work will then be uncertain 
by ±50 % due to the uncertain distance. 

4. BURSTS DETECTED 

Our RXTE coverage detected a total of 6 bursts, in- 
cluding two weak ones, listed in Table El One of these, 
MB3b, was observed in a pair where the second is weaker 
than the first, and the other, XB1, could also have been 
preceded by an unobserved normal burst as it occurs at 
the beginning of a time-series. Four of these bursts, in- 
cluding the pair, were observed by Gemini-S and one of 
these also by HST, giving X-ray, UV, and optical cover- 
age of the same burst. The data quality in the UV, and 
especially the optical, is superb, making these the best 
observed optical bursts in any source. The bursts with 
multiwavelength coverage are shown in Fig. 

One of the simultaneously observed bursts, MB3a, 
shows dips in the otherwise smooth decay. This burst 
occurred at orbital phase 0.478, so does not lie within 
the classical dipping phase-range but our X-ray data do 
appear to show additional dipping near phase 0.5 (Pa- 
per III), so it might be due to transient absorption by 





Elapsed time (seconds) 

Fig. 3. — X-ray bursts with simultaneous multiwavelength cover- 
age. These have been arbitrarily shifted in time to fit on common 
axes. Dashed lines indicate non-contiguous coverage. 



intervening material. In support of this, the dip does ap- 
pear most pronounced in the 2-5 keV energy band, and 
is absent in the optical data. Both these characteristics 
suggest that the dips are due to absorption of the direct 
X-rays along our line-of-sight, rather than variations in 
the intrinsic burst luminosity. 

5. BURST ANALYSIS 

5.1. Cross-correlation Functions 

Cross-correlation functions (CCFs) provide a relatively 
simple technique for examining lags between lightcurves, 
and one which is widely used allowing comparison with 
other datasets. We therefore begin with this approach, 
before developing other techniques customized to the 
reprocessed burst problem. The datasets we have are 
almost exactly evenly sampled, but the Gemini-S pho- 
tometry has a somewhat longer sampling interval than 
the RXTE or HST data. Interpolation cross-correlation 
funct ions ijGaskell fc Peterson! [19871 iWhite fc Peterson! 
1994) are thus ideally suited to these data. We show 
all of the CCFs in Fig. 0] The single, strongly corre- 
lated signal makes identification of the peak very easy, 
and its location is well defined. A lag of a few seconds 
is clearly seen, with individual values tabulated in Ta- 
ble |3| All lags are of order a few seconds, as expected 
for light travel times within the binary. All full strength 
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TABLE 3 
Lags measured from CCF. 



Burst ID Wavelengths CCF Lag Phase 



MB1 


X-ray, opt. 


2.92 


0.20 


MB2 


X-ray, opt. 


4.14 


0.38 


MB2 


X-ray, UV 


4.01 


0.38 


MB3a 


X-ray, opt. 


4.09 


0.48 


MB3b 


X-ray, opt. 


4.30 


0.51 



bursts show rather similar structure, and asymmetry be- 
tween the rising and decaying portions. For MB3a and 
MB3b, the lags derived are consistent, but the shapes 
of the CCF are not. The CCF lags do suggest a posi- 
tive correlation between orbital phase and mean lag, but 
the evidence is not compelling (and would vanish if MB1 
were removed from the sample). Nonetheless, this is as 
expected if there is a significant contribution to the re- 
processed signal from the surface of the companion star, 
as the lag from this will be minimum at phase 0.0 and 
maximum at phas e 0.5. This behavior is also cons istent 
with the finding of Sch oembs fe Zoeschingerl l|1990|) that 
the rise time of optical bursts (without simultaneous X- 
ray coverage) was correlated with the orbital phase in 
the same sense. 

Obviously before extensive interpretation, more so- 
phisticated analysis is warranted to verify these results. 
The quality of the CCFs are superb and the major con- 
cern relates to the meaningfulness of the lag derived. 
Fundamental to the CCF approach is the assumption 
that the optical/UV is simply a lagged version of the 
X-rays, but this is clearly not the case in the burst 
lightcurves, and is not expected theoretically. Several 
factors may be important. One is light travel times. 
These will not only produce a mean lag, but will tend 
to smear out the signal, as dif ferent lags are prod uced 
by different reprocessing sites ijO'Brien et al. 1 12002]) . A 
second factor having a similar effect is that there may 
be a finite diffusion time associated with reprocessing; 
if X-rays deposit energy at a significant optical depth 
it cannot be re-radiated instantly, but there will be fur- 
ther lagging and smearing as the energy diffuses outward. 
Finally, the tails of the bursts are clearly different, indi- 
cating non-linear relationships between the lightcurves. 
This occurs because the bursts have a large amplitude 
and substantially change the temperature of both the 
neutron star and the reprocessing site. As a site cools, 
the peak of the emitted spectrum moves to lower ener- 
gies. In the case of the reprocessed light, it moves out 
of the far-UV, and into the optical, having the effect of 
accelerating the far-UV decay rate (relative to the bolo- 
metric decay) and suppressing the optical decay. None 
of these factors are accounted for in cross-correlations, so 
a more sophisticated model is needed. For example, lags 
can be over or under-estimated if the timing c haracteris- 
tics of the two lightcurves differ l)Koenll2003f) . reflecting 
the more general problem that in this case there is not 
a uniquely defined lag. There is a range, from which a 
CCF favors an average value. 

5.2. Maximum Entropy Deconvolutions 



We next attempt to deconvolve the lightcurves to ob- 
tain the transfer function between X-ray and optical 
wavelengths independent of the width of the burst. This 
will be more readily interpreted than the CCF which is 
much broader than the transfer function. We initially 
use the maximum entropy method (MEM; IHorne et alJ 
I199H lHorndir994() . which does not require an assumed 
model for the transfer function, beyond an assumed de- 
fault used in calculating entropy. In principle, the MEM 
method can then resolve distinct sub-structures, if they 
are present, for example the disk and companion star. 

For each burst we use our X-ray flux lightcurves as the 
driver signal and attempt to fit the optical or UV echo 
over the range —50 to +300 s relative to the burst rise. 
The transfer function was calculated for lags of —50 s 
to +200 s. We use a narrow Gaussian default to min- 
imize enforced smoothing of the derived transfer func- 
tions. The results are shown in Fig. [5J Negligible re- 
sponse was seen before —10 s, and slow declines after 
+50 s; these regions are not shown. 

As expected, the derived transfer functions are much 
narrower than the CCFs. All bursts suggest a single- 
peaked response peaking at a lag of a few seconds, very 
close to the locations of the maxima of the correspond- 
ing CCFs. Besides the main peak, all the transfer func- 
tions have a tail extending to larger lags. This may be 
a real effect, as models of reprocessing of bursts by stel- 
lar atmospheres do predict a small amount of the energy 
emer ges with large delays (jCominskv. London, fc Kleinl 
Il987t see Section T7.2I for more discussion). We cannot 
say with confidence that this effect is real, however, as 
there is another explanation for it. As noted earlier, the 
bandpass-limited response is non-linear, and in particu- 
lar the optical burst decay is much slower than that in 
X-rays because the reprocessed spectrum shifts to longer 
wavelengths (into the optical band) as it cools. A method 
which assumes a linear response will require an extended 
tail to reproduce this. In support of this interpretation, 
the effect is much weaker in the transfer function derived 
from UV data as expected. 

This effect complicates the interpretation of the shape 
of the transfer functions using this method. In addition, 
while it is tempting to ascribe the narrower transfer func- 
tion from MB1 to a phase-dependent effect, the response 
to MB3a is also narrower, and in general one has to be 
cautious in interpreting widths derived from MEM echo 
maps as a broader, smoother peak inherently has higher 
entropy. We therefore seek a less biased approach. 

5.3. Gaussian Transfer Function Fitting 

5.3.1. Method 

iHvnes et all ()1998(1 demonstrated a method for con- 
straining the transfer function by parameterizing it as a 
Gaussian function. Our MEM reconstructions indicate 
that the data do not require multiple resolved compo- 
nents in the transfer function. This does not mean that 
the true transfer function does not include multiple com- 
ponents; only that they are not well constrained by the 
data. This is probably a consequence of the relatively 
long burst timescale, which limits our sensitivity to fine 
temporal structure in the response. Consequently the 
single Gaussian transfer function is an adequate approx- 
imation to this problem and it lends itself well to adapt- 




Fig. 4. — Cross-correlation functions for the bursts. The optical or UV broad-band data are cross-correlated against the integrated X-ray 
flux lightcurves. For MB2, the solid line shows the result with optical data, the dotted one that with UV. For MB3, the solid lines shows 
MB3a, the dotted one MB3b. 




Fig. 5. — MEM reconstructed transfer functions for the simultaneous bursts. For MB2, the solid line shows the result with optical data, 
the dotted one that with UV. For MB3, the solid lines shows MB3a, the dotted one MB3b. 



ing to the rather different assumptions needed for burst 
mapping. 

These arise because of the large amplitude of the vari- 
ations. We have already discussed one problem this 
causes; the X-ray lightcurves are strongly energy depen- 
dent, so determining the driving lightcurve is non-trivial. 
This arises because the X-ray flux arises from heating and 
cooling of the neutron star surface; as the surface cools, 
the spectrum softens and moves out of the higher energy 
bands. We have addressed this by calculating a 2-20 keV 
integrated flux lightcurve, using the spectral information 
to ensure the correct internal relative weighting. Similar 
difficulties exist with the optical and UV data. The opti- 
cal burst is more prolonged than the UV one, for the same 
reason as the bursts are more prolonged at low energies; 
the reprocessed spectrum is moving out of the UV band 
into the optical as the reprocessing regions cool. The net 



effect of this is that the optical and UV fluxes do not re- 
spond linearly to the X-ray flux, and the amplitude of the 
variations is too large for them to be treated as a small, 
linearized perturbation. This is actually a blessing in dis- 
guise however, as we can then use the relative changes 
between the optical and UV fluxes to constrain the tem- 
perature evolution, and hence also the emitting area of 
the reprocessing region. Without this information, there 
would be a degeneracy between temperature and area, 
and indeed this is a problem when we t ry to fit optical 
or UV data independently. For example. [Lawrence et alJ 
(1983) attempted a similar analysis of 4U 1636-536 us- 
ing U BV photometry and were able to estimate that the 
bursts involved a rise in temperature of the reprocessor 
from - 25, 000 K to - 50, 000 K. 

We formulate the problem as follows. We assume that 
an input X-ray lightcurve (the 2-20 keV integrated flux) 
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is a suitable proxy for the bolometric illuminating lu- 
minosity, Lx- Reprocessing is assumed to occur within 
a fixed area A (assumed to be of negligible geometric 
depth) with uniform temperature T(t). The spatial po- 
sition and extent of this region introduces a mean lag, 
t, and a blur a T in the reprocessed bolometric luminos- 
ity relative to the illuminating luminosity. The temper- 
ature is then allowed to vary (with fixed reprocessing 
area) such that T 4 (t) oc L Icp , where the reprocessed lu- 
minosity is assumed to be related to the X-ray luminos- 
ity via the transfer function, , 5: L rcp oc "J * Lx- Given 
the temperature evolution and area we can then pre- 
dict the bandpass-limited UV and optical flux evolution. 
This requires a model for the reprocessed spectrum, for 
which we assume a black body in the absence of a better 
choice. The problem can be parameterized in terms of r 
and a T , defining the transfer function (which is arbitrar- 
ily normalized), T m j n and T m which are the minimum 
temperature and the peak irradiation temperature, de- 
fined such that T^ ax = T^ in + Tj 4 ., and A, the projected 
cross-sectional area of the reprocessing region. Ideally 
this gives enough parameters to fit the UV and optical 
lightcurves simultaneously. In practice, some additional 
nuisance parameters are needed. Since the HST abso- 
lute timing is uncertain, we fit the optical and UV lags 
independently; the optical should be reliable, but the UV 
one is not, comprising both a true lag and a clock uncer- 
tainty. We also allow the optical and UV bursts to have 
different reprocessing areas. This is initially assumed to 
reflect uncertainties in the calibration of the optical pho- 
tometry and in dereddening of the opticakUV flux ratio, 
rather than physical differences. In the ideal case, both 
lags and both areas would be equal. 

5.3.2. The multiwavelength burst (MB2) 

We begin with the burst which has both UV and op- 
tical coverage as this should be best constrained. It also 
allows us to test the reliability of the optical fit by com- 
paring it with the joint optical-UV fit. 

We search for the best fitting m odels (in the \ 2 se nse) 
using a downhill simplex method i|Press et al]ll992|) . To 
guard against converging on a local minimum each fit 
used three passes, each starting at the previous best fit, 
and we also varied the initial starting simplex to ap- 
proach the minimum from different directions in param- 
eter space. We find that the method can reproduce both 
the optical and UV lightcurves very well. The best fitting 
joint model (MB2-1) is shown in Fig.Eland Tabled with 
approximately the same parameters found for a variety 
of starting simplexes. There was a small jitter depen- 
dent on the starting point, but this was much less than 
the uncertainties that we estimate, so all fits agreed to 
within errors. The best fitting reduced x 2 is 1.85, but 
this is pessimistic. The reason is that we have used the 
X-ray lightcurve to construct the model, and its errors 
have not been accounted for. To compensate for this in 
deriving estimates of errors on the fit parameters, we use 
the common technique of rescaling the errors so that the 
best fit model has x 2 /dof = 1. 

The model derived has reasonable parameters. The 
heated region has a quiescent temperature of 18,500 K 
and is heated to a peak temperature of 36,000 K. The 
mean lag is 4.0 s, and there is considerable smearing, with 
a FWHM of 5.8 s. The implied error in the HST timing 
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Fig. 6. — Single black body fits to MB2, allowing the optical- 
to-UV areas to float independently (dashed line), or fixing them 
to the same value (dotted line). The underlying non- burst X-ray 
luminosity that has been subtracted off is 2.1 X 10 36 ergs -1 . 



is only 0.2 s, which is better than might be expected. 
The lag is of the order expected for a combination of 
disk and companion star reprocessing, but the smearing 
involved is large. The inferred emitting projected area 
is 0.9 x 10 21 cm 2 from the UV data and 1.9 x 10 21 cm 2 
using the optical. The areas derived are fully consis- 
tent with expectations from reprocessing on the accre- 
tion disk (1.5 ± 0.2 x 10 21 cm 2 ) and/or companion star 
(2.0 ± 0.9 x 10 21 cm 2 ). The large discrepancy in the area 
derived from the optical and UV is more problematic, 
however, indicating that the optical/UV flux ratio is very 
different from that implied by the temperatures; the op- 
tical is a factor of ~ 2 brighter than expected. This 
is much larger than can be accommodated from errors 
in the dereddening (16%), UV flux calibration (4%) or 
optical calibration (probably less than 10%). The dif- 
ference is clearly significant; a similar fit with the areas 
forced to be equal (MB2-2) results in a visibly poor fit 
to the UV lightcurve. Since the derived magnitudes are 
already at the low end of the range observed by previous 
studies it is more likely that the optical flux has been 
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TABLE 4 

Fits to bursts with Gaussian transfer functions. 



Burst model r opt (s) r uv (s) a T (s) T min (K) T irr (K) A opt a (cm 2 ) A uv a (cm 2 ) X 2 /dof 

MB1-1 2.57 ±0.02 - 0.43 ± 0.01 [18500] 34800 ± 500 (1.88 ± 0.03) x 10 21 - 1.19 

MB1-2 2.44 ±0.04 - 1.23t°;^ [18500] 35400 ± 400 (1.72 ± 0.02) x 10 21 - 1.33 

MB2-1 4.05 ±0.07 3.80 ±0.18 2.48 ±0.08 18500 ±700 35300 ± 1200 (1.92±q$) x 10 21 (0.87±o;Jg) x 10 21 1.85 

MB2-2 4.12±0.08 3.81 ± 0.16 2.62 ± 0.09 13600 ± 100 b 26700 ± 200 b (2.65 ± 0.05) x 10 21 - 2.12 

MB2-3 4.05 ± 0.07 3.80 ± 0.17 2.48 ± 0.08 12000 ± 700 23500 ± 1200 (2.32 ± 0.01) x 10 21 (1.3lj;g;°§) x !0 21 1-85 

MB2-4 4.08 ± 0.07 3.81 ± 0.16 2.54 ± 0.09 9100 ± 100 b 18700 ± 100 b (3.03 ± 0.02) x 10 21 (3.03 ± 0.02) x 10 21 1.90 

MB2-5 3.78 ±0.15 2.43 ±0.21 [18500] 35900 ± 1500 - (0.82±g;J|) x 10 21 1.11 

MB2-6 4.05 ± 0.08 - 2.48 ± 0.10 [18500] 35300 ± 400 (1.92 ± 0.05) x 10 21 - 2.69 

MB3-1 4.02 ±0.09 - 2.69 ± 0.11 [18500] 32700 ± 500 (2.50l°;°4) x 10 21 4.57 



a The errors quoted on areas only include statistical uncertainties; 
distance, reddening, and calibration uncertainties will dominate. 

b These uncertainties are unrealistically small. In this model, the 
temperature is set precisely by the ratio of optical to UV flux, but 
the uncertainty only includes statistical errors. 



underestimated than that it has been overestimated by a 
magnitude as required. If the extinction curve is far from 
the Galactic average, then a larger reddening correction 
to the UV is possible, and this could rectify the problem. 

Alternatively, there may be a deficiency in the model. 
The most likely interpretation is that the reprocessing 
region is not isothermal; indeed this is not expected to 
be the case. The disk is expected to span wide range 
of temperatures from rim to center, and the compan- 
ion star will provide an additional reprocessing site. A 
multi-temperature spectrum will be broader than a single 
temperature one, and if the peak is in the UV (as implied 
by the derived temperatures), then this will flatten the 
optical tail and increase the optical fluxes. A simple test 
of whether this explanation can help resolve the discrep- 
anc y is to use simple i rradiated disk spectra as discussed 
by l|Hvnes et al.ll2002|) . Instead of using black bodies of 
temperature T, we use irradiated disk spectra with irra- 
diation temperature at the edge T, and viscous temper- 
ature set to zero. In practice the latter constraint only 
means that the temperature is dominated by irradiative 
heating; increasing the viscous temperature reduced the 
quality of the fit. Models MB2-3 and MB2-4 show pa- 
rameters derived using this alternate spectral model for 
floating and fixed areas respectively. Clearly this model 
does represent an improvement, as the discrepancy be- 
tween optical and UV areas has been reduced when they 
are allowed to float, and the quality of fit is substantially 
better if they are fixed to be the same. This does not 
appear to be the whole story, however, as a discrepancy 
remains. This is probably a consequence of the cooler re- 
sponse from the companion star (which is inferred to con- 
tribute; see Section We tried adding another black 
body component, but found that beyond this point fits 
were poorly constrained and multiple solutions were pos- 
sible, all yielding good fits with no discrepancy in optical 
and UV areas. We feel that a better way to approach this 
problem i s to use a model of th e binary following the ap- 
proach of O' Brien et al.l l|2002j) . in which the binary ge- 
ometry reduces the independence of the parameters that 
we currently have using arbitrary components. Such a 
model would, however, benefit from information gleaned 
from other aspects of our dataset, for example the orbital 



lightcurves, so we defer this treatment to Paper III. 

As a precursor to examining MB1 and MB3, for which 
only optical data are available, we also tested fitting opti- 
cal and UV lightcurves separately. We initially tried fit- 
ting with all the parameters left free, but found that the 
problem was then very poorly constrained; a single band- 
pass did not allow us to uniquely determine both mini- 
mum and maximum temperatures. The formal minimum 
for an optical-only fit occurs for very low T m j n ~ 3500 K. 
While this may be of order the stellar temperature, the 
regions of the star exposed to bursts are also exposed to 
persistent radiation, and should not be this cool. Fur- 
thermore, this solution is strongly inconsistent with the 
UV lightcurve, and is only a shallow minimum; x 2 is vi r ~ 
tually constant for all < T min < 50, 000 K. We thus 
chose to instead fix T min to the value determined by the 
joint fit, 18,500 K. The parameters derived in this case 
(MB2-5 and MB2-6) were not significantly different from 
those based on the joint fit. The lags cannot be reliably 
compared, due to the uncertainty in HST timing, but 
there is no significant difference between the widths de- 
rived. 

One important statement can be made from these anal- 
yses. Independently of whether we fit optical and UV 
data jointly or separately, and regardless of the spectral 
model assumed, we derive essentially the same lags and 
widths, i.e. the transfer function derived is robust and 
not sensitive to these assumptions. This means we can 
usefully compare transfer functions from different bursts 
that do not have UV coverage, and without knowing the 
correct spectral model to use. 

5.3.3. The double burst (MBS) 

We now consider the double burst, MB3. This is a 
more complex case than MB1 but occurred only about 
four hours (one binary orbit) after MB2 on the same 
night, so there are less likely to have been large changes 
in the accretion flow geometry than for the first burst 
which occurred about 20 hours (5 binary orbits) earlier. 
Since MB3a and MB3b are so close together, we per- 
form a fit to both simultaneously, as this will provide a 
visual check of the consistency between the two bursts. 
The major difficulty with these bursts is that they exhibit 
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dips, most prominently around 35 s and 360 s. These dips 
are clearly not present in the optical data, and appear 
more prominent at low energies (Fig. |3J) . This suggests 
that the dips are due to absorption of the X-rays, and 
should be excluded from any attempt to model the op- 
tical data. We do this by masking out regions of the 
optical lightcurve that would be affected. The masked 
regions are indicated in Fig. 

As for MB2, T m j n is not meaningfully constrained by 
the data in the absence of UV coverage. For T m - ln > 
6000 K, x 2 is virtually constant, and while cooler solu- 
tions are formally favored, we choose to ignore them for 
the reasons described above. A fit with T min = 18, 500 K 
is not visibly worse than one with very low T m j n , so we 
choose to use this value for lack of better information. 
This will facilitate a more direct comparison with MB2. 

The fits to the lightcurves of both MB3a and MB3b 
are shown in Fig. [7J and parameters are given in Ta- 
ble Given that only T m i n is externally constrained, 
the agreement with MB2 is striking. r opt values agree 
within uncertainties and er T virtually agrees. A slightly 
larger width is derived from MB3, but the difference may 
not be significant. It would be expected if the response 
arises from a combination of disk and companion star, 
however, as MB3 occurred around phase 0.5 when the 
lag from the companion star is maximized. The differ- 
ence in temperatures and areas is larger, but may be 
related to each other. We have artificially fixed T m in, 
which will have the effect of also constraining T max . If 
T m i„ had actually increased between bursts, then the cor- 
responding Tj rr would also have been higher, and a lower 
area would have been required to fit the MB3 fluxes. For 
example, a fit with T min fixed to 21,000 K greatly reduces 
the discrepancy between T m and A opt between MB2 and 
MB3. Without UV data, however, we cannot say which 
of these (or other) combinations are correct, so have cho- 
sen to assume that T m i n does not change. 

5.3.4. The first burst (MB1 ) 

We have left MB1 until last as this is the least con- 
strained case. We do not have UV data allowing inde- 
pendent fits of Tmin and T max , but neither are we well 
justified in assuming similarity to MB2 or MB3. MB1 
occurred on the preceding night, and also at phase 0.2 
when our visibility of the system was significantly differ- 
ent. We might then expect any of the parameters to be 
different. Nonetheless, we will begin from the same start- 
ing point as for MB3, i.e. assuming that T m ; n = 18500 K, 
as for MB2, and performing an unconstrained fit on the 
other parameters. This burst is more problematic than 
the preceding ones, as there are actually two plausible 
solutions differing mainly in the width of the transfer 
function. We show in Fig. [5] the dependence of normal- 
ized x 2 on t ne Gaussian width, tr T for all three bursts. 
For MB2 and MB3 there was a single, well-defined min- 
imum, and the location is approximately the same for 
these two bursts. This width is clearly not consistent 
with the data from MB1, and a narrower transfer func- 
tion is required. Formally the best fit occurs for a very 
narrow transfer function: MB1-1 as listed in Tabled x 2 
is a rather complex function of a T is this region, and the 
best fit represents a very narrow minimum (hence the un- 
realistically small uncertainties for <r T in model MB1-1). 
These characteristics suggest that the behavior at small 



o> may reflect random correlations in the noise between 
X-ray and optical lightcurves. A secondary minimum 
is present corresponding to a broader transfer function, 
model MB 1-2. This minimum more closely resembles 
those seen for MB2 and MB3, and we suspect that this 
could be the true solution. In support of this, we show 
both models fitted to the data in Fig. Model MB1- 
1 predicts a steeper optical rise than observed, whereas 
MB1-2 gives approximately the correct rise time. It is not 
obvious in what respect MB 1-1 is preferred, also indicat- 
ing that it could owe more to fitting the noise than to 
fitting the real structure. We will therefore adopt MB 1-2 
as the preferred fit to this burst. 

The parameters obtained do seem plausible. We de- 
rive a shorter lag and narrower lag distribution (inde- 
pendently of whether we adopt MB1-1 or MB1-2). Both 
are to be expected at phase 0.2, as the companion star 
will make less contribution, and will do so at a shorter 
lag. The temperature and area derived are comparable, 
although a smaller reprocessing area is found, again as 
expected at this phase. Based on these parameters there 
is no need to invoke significant differences in the geom- 
etry or temperature of the reprocessing regions between 
the two nights. 

Nonetheless we did attempt fits with a range of T m i n 
values. As for other bursts we find that formally the best 
fits are for unphysically low base temperatures (T m i n < 
4000 K, but that a large range of values larger than this 
also give acceptable fits. There is no reason to prefer a 
different T min value to MB2 or MB3. 

6. THE ULTRAVIOLET BURST SPECTRUM 

As our MB2 UV data were obtained in a spectroscopic 
mode, as well as obtaining the lightcurve of the repro- 
cessed burst, we have a unique opportunity to examine its 
spectrum. Fig. 1101 shows the spectrum of the extra light 
during the first 100 s of the burst. This was constructed 
by extracting a series of 100 s sub-spectra using the STS- 
DAS task INTTAG and processing them in the same way 
as regular spectra. The non-burst spectrum was defined 
from 300 s intervals before and after the burst. More de- 
tails of the spectral reduction will be provided in Paper 
II. 

The burst spectrum appears to be continuum domi- 
nated with only a weak contribution from the lines. C in 
11 75 A is particularly strongly enhanced, and to a lesser 
extent Nv 1240 A, but neither dominate the flux. The 
burst is not pronounced at all in Civ 1550 A or Hell 
1640 A. 

The continuum shape is consistent with the burst mod- 
els described earlier, being well fitted by the difference 
between a black body with temperature ~ 27, 000 K and 
one at ~ 18, 500 K. The upper temperature corresponds 
to the expected flux-weighted average temperature dur- 
ing the 100 s interval contributing to the burst as calcu- 
lated from model 1 for MB2. This shows that the model 
does reproduce the wavelength dependence as well as the 
temporal evolution. It also indirectly supports our low 
reddening value {E{B — V) = 0.06 ± 0.03) derived from 
the 2175 A interstellar feature (Paper III). If this red- 
dening had been underestimated due to an anomalously 
weak 2175 A feature then in general we would not ex- 
pect the temperature implied by the shape of the spec- 
trum to be consistent with that earlier obtained only 
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Fig. 7. — Single black body joint-fit to both MB3a and MB3b. The underlying non-burst X-ray luminosity that has been subtracted off 
is 1.9 X 10 36 ergs~ 1 . Gray bars indicate X-ray dipping regions. These were retained in calculating the model optical lightcurve but were 
masked out of \ 2 calculation. The discrepancy between the model and the optical data at these points is further evidence that the dips 
are a consequence of our changing visibility of the X-ray source, not a true change in the unabsorbed luminosity. 
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Fig. 8. — Normalized x 2 as a function of Gaussian width for each 
burst. For each point the width was fixed, but other parameters 
were optimized to find the best fit. Annotations 1 and 2 for MB1 
refer to the two minima given as solutions in Table [3] Note that 
all curves flatten below <j t ~ 0.3 s, corresponding to a Gaussian 
FWHM less than 0.7 s, and less than the optical time-resolution. 



from time-dependence of the reprocessed light, and the 



far-UV spectrum would be redder. 

7. DISCUSSION 

7.1. Evidence for phase dependence of the burst 
response 

One of the primary goals of echo mapping is echo- 
tomography: phase-dependent echo maps in which dif- 
ferent components, such as the disk and companion star, 
can be disentangled based on different amplitudes and 
phasings of phase dependence. In principle this tech- 
nique could even be used to obtain a direct measure of 
the binary separation and orbital inclination independent 
of other constraints. Real datasets have, however, fallen 
far short of these aspirations, and we certainly would not 
claim that the results presented here fulfill the hopes of 
echo-mapping. Nonetheless, we do see changes in the 
response, and it is worthwhile to investigate if they are 
consistent with expected phase-dependent changes, and 
can constrain the parameters of UY Vol. 

One of the major limitations, of course, is that MB1 
was obtained on a different night to MB2 and MB3. Since 
MB1 provides the main sensitivity to phase-dependent 
changes, we must apply the caveat that the differences 
observed could reflect secular changes in the actual repro- 
cessing geometry rather than just differences in our view- 
ing angle. The persistent (pre-burst) X-ray flux is essen- 
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Fig. 9.— Single black body fits to MB1. The underly- 
ing non-burst X-ray luminosity that has been subtracted off is 
1.6 X 10 36 ergs -1 . The dotted line corresponds to the formal best 
fit, MB1-1 in Table HI The dashed line corresponds to the broad 
secondary maximum which we suspect may be the true solution, 
MB1-2. Note that the latter fits the rise time of the optical burst 
better. 



tially the same for MB1 and MB2. It is lower for MB3, 
but there is some evidence for dipping at this time. This 
implies no substantial difference in the accretion rate be- 
tween the MB1 and MB2 epochs. Optical lightcurves will 
be compared in Paper III and should be very sensitive 
to changes in the reprocessing geometry between epochs. 
The optical lightcurve at the time of MB1 is somewhat 
different, but the differences arc not dramatic and may 
primarily reflect different realizations of the strong flick- 
ering which is present in the optical data. The ampli- 
tude of the differences between MB1 and MB2 epochs 
K 25 %) is comparable to that between MB2 and MB3 
epochs, and also to the amplitude of individual flickering 
events. There is thus no compelling evidence for a change 
in the state of the system, but this possibility cannot be 
securely discounted. 

Assuming it is reasonable to compare the responses 
from the three bursts, the major difference is that MB2 
and MB3 both have responses that are lagged and 
smeared by more. This is as expected from their phases. 
Fig. 1111 shows predicted responses as a function of or- 
bital phase ca l culate d using the methods described by 
lO'Brien et all ipOOl for model 2 (q = 0.2). Given the 
limitations inherent in approximating the true response 
with a Gaussian, the range of lags observed is in rea- 
sonable agreement with predictions. In particular, we 
do expect that MB2 (phase 0.38) will have a similar 
lag distribution to MB3 (phase ~ 0.5), whereas we ex- 
pect a quicker and less smeared response from MB1, as 
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Fig. 10. — Far-UV spectrum of MB2. The upper panel shows the 
spectrum of the first 100 s of the burst compared to that before and 
after. The lower panel shows the difference between the two. The 
dashed line shown in the lower panel is the expected black body 
burst spectrum based on the time-averaged temperature evolution 
fitted earlier. All spectra shown here have been rebinned to 4.7 A 
per pixel. 



observed. The amplitude of the difference between re- 
sponses is in approximate agreement with expectations. 
When plotted in this way, it is apparent that all bursts 
appear to show the onset of a strong response at approx- 
imately the same time, but that MB2 and MB3 extend 
for longer beyond that. 

While the general picture is in agreement with predic- 
tions, it is clear that the data are not precise enough to 
constrain the system parameters further. Nonetheless, 
the very fact that the response appears to change with 
phase, and with such a large amplitude, does indicate 
significant response from the companion, and the large 
smearing of the response around phase 0.5 indicates that 
the disk must also contribute. 

We can also consider these results in the context of 
other observations of this and other sources. The com- 
pariso n of most interest is with lSchoembs fc Zoeschingerl 
( 1990). These authors found that optical bursts (with no 
simultaneous X-ray coverage) exhibited longer rise times 
when the companion star would be on the far side of the 
disk. This is equivalent to saying that the smearing in- 
creases close to phase 0.5 (as we also find), provided that 
the intrinsic rise time of the X-ray bursts is the same. 
The large dispersion in the response that we see close to 
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Fig. 11. — Predi cted echo tomogram calculated according to 
lU'Brien et"aT] 120021) . The broad non-varying region to the left 
is the disk, the moving component is the irradiated face of the 
companion star. Solid bars indicate the response inferred from 
Gaussian fitting of MB1 (lower), MB2 (middle), and MB3a and 
b combined (upper). The bar represents the lags over which the 
Gaussian model is above half response (i.e. its length is the full- 
width at half maximum). 



phase 0.5 is als o not unprecede n ted. Similar dispersions 
were f ound bv iTruemper et all lfl985]l for 4U 1636-536, 
and bv lKong et alJ l|2000ft forGS 1826-24. In all of these 
cases they likely indicate contribution from both disk and 
companion, as a response from the companion alone will 
vary in lag, but should always have a relatively narrow 
response. 

7.2. Local diffusion time delays 

There are two likely sources of smearing of the re- 
processed signal. The delays most commonly discussed 
(and considered so far) are those from light-travel times 
across the binary; global delays. A finite reprocessing 
time could also introduce a local delay. X-rays above 
about 1 keV will experience a low-enough photo-electric 
absorption opacity to deposit energy at a significant op- 
tical depth in an atmosphere. Reprocessed energy will 
then have to diffuse outward, and so there will be an 
additional diffusion time delay. 

Thi s topic was considered initially by iPedersen et alJ 
( 1982) in analyzing bursts from 4U 1636-536. They esti- 
mated that the diffusion time for a typical photon would 
be ~ 0.6 s, largely negligibl e compared to the expected 
light travel time delays. iCominskv. London, fc Kleinl 
( 1987) examined the problem more rigorously, calculat- 
ing time-dependent responses of a hot stellar atmosphere 
to an X-ray burst, including the effect of the burst on the 
atmospheric temperature structure and opacities. Their 
result s were in agreement with those of IPedersen et aTl 
( 1982), and they found that 50 % of the reprocessed light 
is expected within just 0.2 s, but that there was also a 
very extended tail to the response up to 10 s. Finally, 
iMcGowan et alJ l)2003|) also considered this issue, but ne- 
glected photo-absorption, which will be critical for the 
soft X-rays produced in a burst. Based on these cal- 
culations, we would not expect our method to be sig- 
nificantly affected by reprocessing times, as a Gaussian 



response approximation will be dominated by when most 
of the light emerges, and this occurs within a fraction of a 
second. Indeed, our observations are broadly consistent 
with expectations from light travel time delays and, al- 
lowing for the limitations of the Gaussian approximation, 
do not obviously require additional local delays. 

This conclusion may be different from that obtained 
for correlated flickering in the l ow-mass X-ray b i naries 
Sco X-l and LMC X-2, for which IMcGowan et all pOOl 
argue that lags may be too long for light travel times 
alone. More and better data are needed to confirm 
this discrepancy, but if confirmed the key difference may 
be i n the nature and spectrum of the X-ray irradia- 
tion. ICominskv. London, fc Kleinl l)1987|) also considered 
harder irradiation than that provided by an X-ray burst, 
and found that it could increase the di ff usion timescale 
significantly (see also IMcGowan et~aT] pOO^) ). Thus 
while reprocessing times are of marginal significance in 
considering burst reprocessing, they may be of more im- 
portance for other applications where further investiga- 
tion is needed. The corollary of this is that the relatively 
soft irradiation of bursts may provide the least biased 
way to isolate light travel time delays and hence the ideal 
signal for echo-tomography. 

7.3. Expected temperatures 

It is interesting to compare the deduced reprocessing 
temperatures with those expected based on the X-ray 
luminosity, as this provides a test of the reprocessing 
efficiency. We will consider MB2, for which the tempera- 
tures are best constrained. For this burst, the persistent 
luminosity (at 6.8 kpc) was 2.1 x 10 36 ergs -1 and the peak 
burst luminosity was 6.9 x 10 37 ergs — x . 

Irradiation of the companion star is sensitive to the 
mass ratio, as smaller companions will more readily be 
shielded by the disk. For example, if we take model 1 
(q = 0.08) and th e typical effective disk opening angle 
of 12° deduced bv Ide Jong, van Paradiis. fc Augusteiinl 
(1996) then the companion would be completely shielded. 
This is clearly not the the companion can di- 

rectly eclipse the neutron star, but we cannot quantify 
how directly it is illuminated. Note that the opening 
angle need not be that of the disk rim, but could rep- 
resent material in the inner disk that can shield the 
companion. For models 2 and 3, calculations are more 
straightforward and the differences between the two are 
less dramatic. For a disk opening angle of 6° we ex- 
pect angles of incidence of greater than 45° from nor- 
mal, increasing to 60° for opening angles of 12°. We 
assum e an X-ray albedo for the companion of ~ .4 fol- 
lowing dc .long , van Paradiis. fc Augusteiinl l)199fitL We 
then expect the highest persistent irradiation tempera- 
ture on the companion of 16, 000-21, 000K depending on 
mass ratio (for the range 0.20-0.34) and disk opening 
angle (in the range 6-12°). This is the temperature at 
the point closest to the compact object and most of the 
irradiated region is at temperatures less than this as the 
angle of incidence is steeper. At the burst peak we would 
expect temperatures of 38,000-50,000 K. 

For the disk we estimate the estimated temperature 
at the disk r i m usi ng th e semi-empirical pr escription of 
IDubus et all l|1999l) and iDubus et alJ l|2001fl . Assuming 
Dubus' estimate of the irradiation efficiency and a disk 
radius of 5 x 10 10 cm we expect a persistent irradiation 
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temperature of 9000 K and a peak burst temperature of 
21,000 K. Both of these are the lowest temperatures in 
the disk, and most of the disk area is hotter than this 
(the opposite case to that calculated for the compan- 
ion star). For comparison, the persistent luminosity of 
2.1 x 10 36 ergs _1 corresponds to a mass transfer rate of 
~ 1.2 x 10 16 gs _1 . For a steady state lobe-filling disk the 
effective temperature due to viscous heating should be 
just 2500 K, indicating that the disk is in a regime that 
can only remain in a high state with irradiative heating. 

Our earlier conclusions that the reprocessed bursts 
likely arise from a combination of emission from the com- 
panion star and disk are consistent with these calcula- 
tions. Our estimate of the peak reprocessing tempera- 
ture is 35,000 K. The irradiated area of the companion 
star is expected to have temperatures possibly extend- 
ing up to 50,000 K, while the disk could plausibly have 
temperatures from 21,000 K upward. The temperature 
ranges are overlapping, with the greater distance to the 
companion offset by its more direct illumination, and are 
consistent with our observations. 

8. CONCLUSIONS 

We have reported several simultaneous X-ray bursts in 
the low-mass X-ray binary UY Vol, including the high- 
est quality reprocessed optical coverage, and the only re- 
processed UV burst that we are aware of in any source. 
These results allow a more thorough test of the paradigm 
for reprocessing of X-ray bursts, and reprocessing of X- 
ray emission in general, than previously possible. Several 
key assertions about reprocessing have been tested. 

i) X-rays are absorbed at relatively high optical depths, 
thermalized, and re-emitted with a quasi-black body 
spectrum. By obtaining the UV spectrum of the ex- 
tra light produced during a burst, we have shown that it 
is indeed dominated by continuum emission, and that 
the shape of that continuum is consistent with black 
body emission with temperatures as inferred from the 
lightcurves of the bursts. The discrepancy with the op- 
tical flux, however, suggests that a single-temperature 
black body model is not sufficient. 

ii) Reprocessed emission in the optical and UV is 
lagged predominantly by light travel times rather than 
by local photon diffusion. This appears to be borne out 
by our observations as light travel times are sufficient to 



Chabrier, C, & Baraffe, I. 2000, ARA&A, 38, 337 
Cominsky, L. R., London, R. A., & Klein, R. I. 1987, ApJ, 315, 
162 

de Jong, J. A., van Paradijs, J., & Augusteijn, T. 1996, A&A, 314, 
484 

Dubus, C, Lasota, J. -P., Hameury, J.-M., & Charles, P. 1999, 

MNRAS, 303, 139 
Dubus, G., Hameury, J.-M., & Lasota, J.-P. 2001, A&A, 373, 251 
Fitzpatrick, E. L. 1999, PASP, 111, 63 

Fukugita, M., Shimasaku, K., & Ichikawa, T. 1995, PASP, 107, 945 
Gaskcll, C. M., Peterson, B. M. 1987, ApJS, 65, 1 
Gottwald, M., Haberl, F., Parmar, A. N., & White, N. E. 1986, 
ApJ, 308, 213 

Haswcll, C. A., Hynes, R. I., King, A. R., & Schenker, K. 2002, 

MNRAS, 332, 928 
Home, K. 1994, Astronomical Society of the Pacific Conference 

Series, 69, 23 

Home, K., Welsh, W. F., & Peterson, B. M. 1991, ApJ, 367, L5 
Hynes, R. I., O'Brien, K., Home, K., Chen, W., & Haswell, C. A. 
1998, MNRAS, 299, L37 



explain the lags and smearing observed. Note that this 
conclusion (and the preceding one) may not be true in 
general, but may be specific to irradiation with the spec- 
trum of an X-ray burst. Harder or softer irradiation is 
likely to produce a different response. 

iii) X-rays are reprocessed both by the accretion disk 
and by the companion star. This is supported by our 
observations, as we do appear to see a phase-dependence 
of both the lag and the smearing of the response. If the 
response came from the disk alone, we would expect little 
variation. If it came from the companion alone, we pre- 
dict variations in the lag, but not in the smearing. The 
temperatures that both the disk and companion are ex- 
pected to reach during a burst are consistent with our ob- 
servations. The failure of a single-zone reprocessor model 
to explain the relative UV and optical fluxes in the burst 
also supports a multi-component response. Our observa- 
tions do not appear to resolve distinct lags from disk and 
companion, however, possibly a consequence of the long 
bursts exhibited by UY Vol. Sources which exhibit short 
duration bursts may be better suited to this analysis. 
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